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It is studied the bright spatial solitons in nonlocal defocusing Kerr media with parity-time (PT) symmetric 
potentials. We find that these solitons can exist and be stable over a different range of potential parameters. 
The infiuence of the degree of nonlocality on the solitons and the transverse energy fiow within the stable 
solitons are also examined. © 2011 Optical Society of America 
OCIS codes: 190.3270,190.6135. 



Generally, defocusing Kerr nonlinearity can induce an 
enhanced beam broadening so that it does not support 
any localized structures other than vortex and dark soli- 
tons, which require background beams [1,2]. However, 
the bright solitons, which were thought to exist only for 
focusing nonlinearity [1], can in fact exist in periodic 
photonic lattices or waveguide arrays with defocusing 
nonlinearity [3-5] . This means that the strong transverse 
periodic refractive index potentials can not only suppress 
the beam broadening due to diffraction but also over- 
come the broadening effect due to the defocusing Kerr 
nonlinearity [1]. 

The nonlocal nonlinear behavior in a system is related 
to the change of nonlinear refractive index, which de- 
pends not only on the intensity of a local beam but also 
on the intensity of the surrounding region due to a spe- 
cific distribution [6] . The nature and extent of nonlocal- 
ity substantially depend on the materials. Principally, 
new effects attributed to nonlocality have been studied 
in thermo-optic media [7], photorefractives [8], and liq- 
uid crystals [9, 10]. 

The definition of PT potentials and their proper- 
ties were discussed in the past few years [11-13]. The 
real part of PT symmetric potentials must be a sym- 
metric function of position, while the imaginary com- 
ponent should be antisymmetric. Recently, parity- time 
symmetric potentials have been introduced into optical 
field [14-27]. However, thus far all studies focus on bright 
solitons in self-focusing optical PT symmetric media, and 
bright spatial solitons in nonlocal Kerr self-defocusing 
media with a single PT complex potential are never re- 
ported. 

In this paper, we investigate the bright spatial soli- 
tons in nonlocal defocusing Kerr media with parity-time 
(PT) symmetric potentials. It is found that these self- 
tapped states can exist and be stable over a different 
range of potential parameters. In addition, we show the 
influence of the degree of nonlocality on the solitons and 
the transverse energy flow within the stable solitons. 



In a nonlocal Kerr self-defocusing medium with PT 
symmetric potentials, the one-dimensional optical beam 
evolution is governed by the following normalized non- 
linear Schrodinger-like equation for q and 0, which are 
respectively the dimensionless light field amplitude and 
the nonlinear correction to the refractive index [14-25], 



dz 



^[V{x)^iW{x)]q-H = ^. 
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where z is the propagation distance, V{x) and W{x) are 
the real and the imaginary components of the complex 
PT symmetric potential, respectively. V{x) is an even 
function and W{x) is odd. Physically, V{x) is associated 
with index guiding while W{x) represents the gain/loss 
distribution of the optical potential. We assume that the 
depth of the refractive index modulation is small com- 
pared to the unperturbed index. The nonlocality of the 
materials is supposed to be ruled with an exponential re- 
sponse function g{x) = 1 / {2a^ ^ eyi'^{—\x\ / ^ (as in 
liquid crystals), where a is the degree of the nonlocality. 
We are going to search for a stationary soliton solution 
of Eq. (1) in the form of q{x^z) = u{x)e^^^ ^ where u is 
a complex function and h is the propagation constant of 
spatial solitons. In this case u satisfies 
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dx'^ 

Here, we assume a Scarff II potential shown in 
Fig. 1(a) where V{x) = V^sechixY and W{x) = 
Wosec/i(x)tanh(x), with Vo and I^o being the amplitudes 
of the real and imaginary part [12, 14, 15]. Although the 
PT symmetric potential has crossed the phase transition 
point, the solitons still exist because the amplitude of the 
refractive index distribution can be altered by the beam 
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itself through the optical nonlocal nonlinearity. The PT 
symmetric will remain broken if it cannot be nonlocal 
nonlinearly restored [14, 15, 17, 18]. 

To check the stability of the solitons with the method 
of linear stability analysis, we assume q{x) = u{x)e'^^^ + 
e[F(x)e*^^ +G*(x)e-*'^*^]e*^^ where e < 1, F and G are 
the perturbation eigenfunctions, and 6 is the growth rate 
of the perturbation. By linearizing Eq. (1), we gain [14] 

SF = + (V + iW)F -^nF-bF^ uAn, (3a) 

5G = -—r + (-V + iW)G -nG^bG- ii*An. (3b) 

where n = — g{x — X)\u{X)\'^dX, and An = 

-f^^g{x - X)[G{X)u{X) + F(A)ii*(A)]dA. The bright 
spatial soltions are linearly unstable when S has an imag- 
inary component, on the contrary, they are stable when 
S is real. 
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Fig. 1. (color online) (a)-(d) The real (solid red curve) 
and imaginary components (dash blue curve) of the soli- 
ton solutions: (a) b = 1.288, a = 2, (b) 6 = 1.481, a = 2, 
(c) b = 1.141, a = 5, and (d) b = 1.599, a = 5. (e) 
Transverse energy flow S of solitons (blue and red curves 
represent the solitons in Figs. 1(b) and (d), respectively.) 
(f) lm{S) versus 6 at a = 2 (solid blue curve) and at 
a = 5 (dash red curve), where the points marked with 
circles correspond to the solitons in Figs. l(a)-(d), and 
the inset depicts the curves from b = 1.45 to b = 1.55. 
The potential parameters are Vq = 3.01 and Wo = 0.3. 

In order to gain the solitons solutions, we numerically 
solve Eq. (2) using spectral renormalization method [28]. 
We find a family of the simplest ground-state nonlocal 
bright solitons with PT symmetric potentials. To illus- 
trate the properties of the solitons, we vary 6, a. Wo, and 
Vq. First, the typical cases of the solitons with different 
b and a are shown in Fig. l(a)-(d), where Vq = 3.01 and 
Wo = 0.3. Here the fields of unstable and stable soliton 
are respectively depicted. Evidently, the change of b and 
a can influence the solutions. To shed more light on the 
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Fig. 2. (color online) (a)-(d) The real (solid red curve) 
and imaginary components (dash blue curve) of the soli- 
ton solutions: (a) Wo = 0.38, Vq = 3.01, (b) Wo = 0.60, 
Vo = 3.01, (c) Vo = 2.8, Wo = 0.3, and (d) Vq = 3.15, 
Wo = 0.3. (e) and (f) Transverse energy flow S of solitons 
in Figs. 1(a) and (d), respectively. The other parameter 
is a = 5. 



properties of the stable solitons, we study the parame- 
ter S = {i/2){uu* — u*Ux)^ which is associated with the 
transverse power flow density or Poynting vector across 
the beam [14]. The transverse energy flow density S is 
shown in Fig. 1(e), and the fact that S is positive ev- 
erywhere implies that the energy always flows in one di- 
rection, i.e., from the gain region toward the loss region 
[14]. It is shown in Fig. 1(f) how Im((5), the imaginary 
component of varies versus b. Since Im((5) is not always 
equal to zero, the solitons will not always be stable. One 
can see that, when a is different, the stable regions of 
the solitons are different too, which means the degree of 
nonlocality can influence the stability of solitons. 

Next, we set a = 5 and then vary Wo and Vq. In 
Figs. 2(a)-(d), the bright solitons in different conditions 
are shown. Figs. 2(e) and (f) depict the transverse energy 
flow S of the solitons in Figs. 1(a) and (d), respectively. 
Similarly, one can find that S is positive everywhere, 
which also means that the energy always flows in one di- 
rection. When Vq = 3.01, we gain the relations of b and 
Im((5) with Wo shown in Figs. 3(a) and (b), respectively. 
Fig. 3(a) shows that the bright solitons can exist when 
< Wo < 3.13. Thus for a fixed value of Vq, there exist 
a threshold for the imaginary amplitude Wo. Above this 
threshold (Wonpp = 3.13), a phase transition occurs and 
the spectrum enters the complex domain [14]. Moreover, 
Figs. 3(c) and (d) show how b and Im((5) vary with Vq 
when Wo = 0.3. In Fig. 3(c), one finds that the bright 
solitons can exist when Vo > 1.03, and thus, for given 
Wo, this effective potential nonlinearity will shift the PT 
Vq threshold {Voiow = 1.03) and in turn allow nonlinear 
eigenmodes with real eigenvalues to exist. In Fig. 3(b) 
and (d), lm{S) is not always equal to zero, so the solitons 
are not always stable. To examine linear stability results. 



2 



1.5 
1 

-Q 0.5 


-0.5 



(a) 











0.5 







-0.5^ 






5 3.05 3.1 


< 





(b) 



Fig. 3. (a) The propagation constant b versus Wq at 
Vo = 3.01. (b) lm{5) versus Wo at Vo = 3.01. (c) The 
propagation constant b versus Vq at Wq = 0.3. (d) Im((5) 
versus at I^o = 0.3. Points marked with blue circles 
correspond to the solitons in Figs. l(a)-(d), respectively. 
Points marked with red circles correspond to the thresh- 
old Woupp = 3.13 for the existence of the bright solitons 
where b becomes negative in (a) and (b). Points marked 
with red circles correspond to the threshold Vqiow = 1-03 
for the existence of the bright solitons where b becomes 
negative in (c) and (d). The insets depict the curves from 
1^0 = 3 to Wo = 3.13 in (a) and (b). The other parame- 
ter is a = 5. 



we further simulate the propagation of beams with ran- 
dom noise in the different conditions corresponding to 
Figs. l(a)-(d) and Figs. 2(a)-(d), see Figs. 4. 

To summary, the bright spatial solitons in nonlocal 
defocusing Kerr media with PT symmetric potentials 
are studied. The existence, stability, and propagation dy- 
namics of such solitons are discussed in detail. 
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Fig. 4. (a)-(d) Simulated propagation of the solitons with 
1% random noise corresponding to Fig. l(a)-(d). (e)-(h) 
Simulated propagation of the solitons with 1% random 
noise corresponding to Fig. 2(a)-(d). 
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